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We measure the intensity correlation function of two interfering spatially displaced copies of a 
phase fluctuating Bose-Einstein Condensate (BEC). It is shown that this corresponds to a mea- 
surement of the phase correlation properties of the initial condensate. Analogous to the method 
used in the stellar interferometer experiment of Hanbury Brown and Twiss, we use spatial intensity 
correlations to determine the phase coherence lengths of elongated BECs. We find good agreement 
with our prediction of the correlation function and confirm the expected coherence length. 

PACS numbers: 03.75.Hh, 03.75.Nt, 39.20.+q 
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Since the first realization of Bose-Einstein condensa- 
tion in dilute atomic gases their coherence properties 
have attracted considerable theoretical and experimen- 
tal interest. This interest is due to the central role of 
the coherence properties for the theoretical description 
and conceptional understanding of BECs and their use 
as a source of coherent matter waves in many promising 
applications. 

Remarkable measurements demonstrated the phase co- 
herence of three-dimensional (3D) condensates well below 
the BEC transition temperature T c 0, 0, Q and even at 
finite temperature 0]. However, low dimensional systems 
show a qualitatively different behavior. In particular, it 
has been predicted that one-dimensional and even very 
elongated, three-dimensional BECs exhibit strong spatial 
and temporal fluctuations of the phase while fluctuations 
in their density distribution are suppressed 0, 0] . Thus, 
the coherence properties are significantly altered, result- 
ing in a reduced coherence length which can be much 
smaller than the condensate length. In this case the de- 
generate sample is called a quasicondensate. This regime 
has been subject of recent theoretical efforts based on the 
Bogoliubov and Popov 00 theory to extend and gen- 
eralize the description of BEC and to calculate correla- 
tion functions [lfj . The phase fluctuations are caused by 
thermal excitations of low energy axial modes and thus 
depend strongly on the temperature and trapping geom- 
etry. In particular, for very elongated BECs a nearly 
phase coherent sample can only be achieved far below 
T c . Phase fluctuations were first observed using the for- 
mation of density modulations during ballistic expansion 
of a condensate 0, 0, El . In addition, their effect on 
the momentum distribution has been demonstrated using 
Bragg spectroscopy [3 0] . Nonequilibrium properties 
of these condensates have been studied using a conden- 
sate focussing technique . 

In this Letter we report on a direct measurement of the 
spatial correlation function of phase fluctuating BECs. 
To measure the phase correlation properties we inter- 
fere two copies of a BEC with a spatial displacement 



d (Fig. [TJ. The measured interference pattern is deter- 
mined by the phase pattern of the original condensate 
and a global phase difference between the two copies 
introduced by the interferometer. By varying the dis- 
placement d the first order correlation function can, in 
principle, be measured. However, this method is very 
sensitive to fluctuations in the global phase difference 
and the measurement is further complicated by the sta- 
tistical nature of phase fluctuations. We show that the 
use of intensity correlations in the interference pattern 
overcomes these problems and provides the desired infor- 
mation about the phase correlations in the initial con- 
densate. The result of this measurement is described by 
the spatial second order correlation function and yields 
the phase coherence length of the BEC. In many respects 
our measurement is closely related to the stellar interfer- 
ometer of Hanbury Brown and Twiss |l7lll8j ]. They mea- 
sured the intensity of starlight falling onto two detectors 
and computed the intensity correlations electronically as 
a function of the detectors' distance. This measurement 
yields the transverse coherence length of the stellar light 
and allowed them to determine the stars' diameter. Un- 
like the Michelson stellar interferometer, which uses first 
order correlations, atmospheric fluctuations do not dis- 
turb this measurement. In the same way, global phase 
fluctuations in the interferometer setup do not disturb 
our measurement. 

Our experiments were performed with 87 Rb Bose- 
Einstein condensates in the \F = l,rnjr = — f) hyperfine 
ground state. Further details of our experimental appa- 
ratus were described previously [l^ . The confining po- 
tential was provided by a cloverleaf type magnetic trap 
with an axial trapping frequency of lo x = 2it x 3.4 Hz and 
a radial frequency adjusted between uj r = 2n x 300 Hz 
and u r — 2n x 380 Hz. The number of condensed atoms 
N was varied between 4 x fO 4 and 6 x fO 5 . To allow 
the system to reach an equilibrium state we typically 
waited for 4 s after obtaining BEC by evaporative cool- 
ing (with rf " shielding" ) (23. The interferometric scheme 
is shown in Fig.^and is based on two ir/2 Bragg diffrac- 
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FIG. 1: (a) The interferometer is realized by two tt/2 Bragg 
diffraction pulses. We apply the first pulse after a ballistic 
expansion time of 2 ms to reduce mean-field effects and atomic 
scattering. The spatial displacement d is determined by the 
time At between the pulses. Typical line density profiles are 
shown for a phase coherence length of L$ ~ 25 /im with d — 
7 /im (b) and d = 35 /im (c). Only one of the two output ports 
is displayed. 



tion pulses. These pulses were produced by two coun- 
terpropagating laser beams with a frequency difference 
set to the two-photon resonance. They were detuned 
by about 3 GHz from the atomic resonance to suppress 
spontaneous emission. The pulse duration of 100 fj,s was 
chosen long enough to avoid higher order diffraction and 
sufficiently short to avoid any sensitivity to the internal 
velocity distribution of the phase fluctuating BECs. Af- 
ter a time-of-flight of 30 to 40 ms the two output ports 
spatially separate and the atoms were detected by res- 
onant absorption imaging. We integrate the absorption 
images along the radial direction of the condensate and 
subtract the thermal background to obtain line density 
profiles. 

Typical measured line density profiles are shown in 
Fig. ^ These profiles provide a qualitative argument 
for the use of intensity correlations as a particulary use- 
ful tool for our measurements. When the displacement 
d between the overlapping clouds is chosen smaller than 
the phase coherence length in the sample (Fig. QJ)) re- 
gions with almost identical phases are brought to overlap 
and the resulting interference signal is rather smooth. If 
however d is larger than the phase coherence length, re- 
gions with substantially different phases overlap, result- 
ing in an irregular but high contrast interference signal 
(Fig- ■ I n both cases an average over many realizations 
and relative phases results in the same smooth profile in 
the two output ports of the interferometer, revealing no 
information about the coherence properties. Nonethe- 
less Fig. clearly contains information about the spa- 
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Numerical computation of the function f^ 2 \s,d) = 
d s+d s=d =s+dj for three different d. From top to 
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bottom: d/L = 0.3, 0.2, 0.1. The particular choice of positions 
corresponds to the experimentally relevant case. 



tial coherence properties. An appropriate analysis of the 
correlations in the density profile yields an intensity cor- 
relation function that does not vanish in an averaging 
process and contains the desired information about the 
coherence properties. 

We therefore start our analysis by calculating the spa- 
tial second order correlation function for the case of p hase 
fluctuating BECs. In its most general definition [H| it is 
given by 



g( 2 \x 1 ,x 2 ,x 3 ,x 4 ) 
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where (• ■ -)t denotes an average over an ensemble at 
thermal equilibrium at temperature T. It contains 
the spatial intensity correlation function g( 2 >(xx, x%) = 
g( 2 \x\, X2, X2, x\) as a special case. For 3D condensates 
with repulsive interactions in elongated trapping poten- 
tials density fluctuations are suppressed by the mean- 
field potential 0, fl3l IHj . Therefore the total field opera- 
tor of the condensed atoms can be written as ip(x) — 
\J uq(x) exp(i(j)(x)) , where 4>{x) is the operator of the 
phase |fjj and tiq(x) is the density in the Thomas- Fermi 
approximation. For this field operator one can show that 
all higher order correlation functions can be expressed as 
products of the first order correlation function [i^J] ■ Us- 
ing the explicit expression of the phase operator given by 
Petrov et al. we obtain the second order correlation 
function: 
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where Pj^ are Jacobi polynomials, j is a positive inte- 
ger, 2L is the condensate length and 



L _ 15N (huj x )' 2 
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(4) 



is the phase coherence length in the condensate center 
(in units of L). Here n denotes the chemical potential 
and ks the Boltzmann constant. The function is 
shown in Fig.|2]and contains the functional form of phase 
fluctuations in elongated condensates. All experimental 
parameters are contained in the phase coherence length 

Let us demonstrate how the interferometric sequence 
described above can be used to measure g^ 2 \ In each 
output port of the interferometer the superposition of 
two ballistically expanded spatially displaced copies of 
the original wave function is produced. The field operator 
of the atoms in one output port can be expressed as 

&(x,d) = \Wn{x - d/2)e^ + (*- d / 2 ) 



(5) 



with ^(x) = 4>{x) ± fix + ax 2 26]. The linear term 
results from the mean-field repulsion between the copies 
of the condensate [2lJ and the quad ratic term originates 
from the self-similar expansion |22l . The self-similar 
expanded density distribution n(x) differs only slightly 
from the initial distribution because the increase of axial 
size during expansion is small for very elongated conden- 
sates (about 1% for our experimental conditions). The 
relative global phase between the two overlapping clouds 
0rci = S12 At+8(f>cS is determined by the Bragg diffraction 
process, where 5<fi e H is the change in the relative phase of 
the Bragg beams between the two pulses and S\2 is the 
detuning from the two-photon resonance. In our exper- 
iment Stfies can be controlled by using an electro-optical 
modulator |2^. 5\2 depends on the frequency difference 
of the Bragg beams and the axial release velocity of the 
condensate. Note that even a small change in the release 
velocity of 0.033 /im/ms leads to a phase chan ge o f 7r/2 
for a time At = 3 ms between the Bragg pulses |27j . The 
insensitivity to such randomly varying global phases is a 
major advantage of the intensity correlation method. 

In analogy to the definition of the correlation coeffi- 
cient we define a normalized intensity correlation func- 
tion 
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with the intensity operator — 'i , ^(xi.2,d)^pf(x\ t 2 1 d). 
Here the averages are taken over an ensemble in thermal 
equilibrium and over all relative phases. The possible 

(21 

values of 7^ range between +1 (perfect correlation) and 



-1 (perfect anticorrelation); if 1\ and I2 are uncorrelated 
7y 2 ' ) = 0. By substituting Eq. (J5J) into Eq. JfjJ) we obtain 

7) (xi,x 2 ,d) — cos[2(xi - x 2 )(ad — (3)\ 

x exp 
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The normalized intensity correlation function is the prod- 
uct of two contributions. A cosine resulting from the self- 
similar expansion and mean-field repulsion and an expo- 
nentially decaying term containing the influence of the 
phase fluctuations. The decay constant of this function 

is the phase coherence length. Comparison with Eq. @ 

(2) 

shows that a measurement of 7^- is equivalent to a mea- 
surement of the second order correlation function of the 
trapped condensate. 

(2) 

The measurement of 7^ is performed in the following 
way: For a given trap configuration, evaporative cooling 
ramp and fixed displacement d, a series of measurements 
is recorded by scanning the relative phase of the Bragg 
beams with the electro-optical modulator in small steps 
between and 2ir. This insures that the global phase 
ro i contains all values with equal probability. We ex- 
perimentally determine 7?-' analogous to Eq. ©. The 

ensemble averages (I(x)) are obtained by averaging all 
interference patterns I(x) recorded in a series. Then 
the quantity I(x) — (I(x)) can be determined for each 
realization. The average of these values according to 
Eq. © yields ■ To simplify the analysis we evalu- 

(2) 

ate 7^ at symmetric positions around the center of the 
interference pattern such that x\ — —x 2 — s/2. Then 

(2) (2) 

7). '(s, d) = 7^ '(—s/2, s/2, d) can be expressed as a func- 
tion of s for a given displacement d. Typical results are 
shown in Fig. EI clearly displaying the functional form of 
a damped cosine. 

To extract quantitative results we fit the measured 
function with the theoretical one given in Eq. 0. The 
fit contains only 1$ and the frequency of the cosine as 
free parameters. Figure^lcompares measured correlation 
functions with the corresponding fits, confirming the ex- 
cellent agreement with the expected functional form. Al- 
though the phase coherence length also depends on the 
trapping potential and the temperature, the data sets 
shown in Fig. [3] differ mainly due to the number of con- 
densed atoms, which was set to 4.4 x I0 4 (»), 2.9 x I0 5 (+) 
and 5.0 x 10 5 (d). Except for the smallest atom number 
a minimum is clearly visible and unambiguously defines 
the frequency of the cosine jl^] . The damping of this os- 
cillation yields the phase coherence length. As expected 
(see Eq. (@J), the damping is stronger for small N . In 
case of the smallest No no oscillation is visible, indicating 
a significant decrease of the phase coherence length. 

We have performed such measurements for a large vari- 
ety of atom numbers and temperatures. The good agree- 
ment between the measured phase coherence length ob- 
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FIG. 3: Measured correlation functions (points) compared 
to fitted theoretical curves (lines). Each measured function 
is based on « 25 realizations and plotted up to s = O.SL. 
Neighboring points are not independent since they are based 
on a common set of experimental realizations. Inset: Nu- 
merical simulation including the phase and density evolution 
(points) and the analytical function Eq. Q (line), both for 
the parameters of the Lj, — 40 )im curve. 



tained from the fit and the theoretically predicted one is 
shown in Fig. By varying the displacement d we have 
confirmed that the measured phase coherence lengths are 
independent of this parameter. In all cases the phase 
coherence length was much smaller than the conden- 
sate length, which ranged from 280 /im up to 420 /im, 
i.e. our measurements were performed in the quasicon- 
densate regime. 

So far we have neglected the evolution of the phase fluc- 
tuations during ballistic expansion. This evolution leads 
to a change of the original phase pattern and the appear- 
ance of density modulations. Using the full evolution of 
the wave function |20j we have calculated the expected 
phase change during time-of- flight to be less than tt/10 
for our parameters. To evaluate the influence of density 
modulations on our measurements we have performed 
numerical simulations including the phase and density 
evolution. The inset in Fig. [21 compares a numerical re- 
sult with Eq. 0. The excellent agreement demonstrates 
that this evolution can be neglected for our intensity cor- 
relation measurements and justifies the use of Eq. J7J) to 
extract the phase correlation properties. 

In conclusion, we have demonstrated a new interfer- 
ometric method which allows us to measure the spatial 
correlation function of phase fluctuating BECs. First the 
second order correlation function was calculated and it 
was shown that this function can be measured using in- 
tensity correlations in the interference pattern. Our mea- 
surements were then compared with these results, con- 
firming both, the expected functional form of the cor- 
relation function and the phase coherence length of the 
sample. We have confirmed that this technique is insen- 
sitive to fluctuations of the relative global phase during 
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FIG. 4: Measured phase coherence length compared to the 
theoretically expected one given by Eq. (JIJ. The dotted line 
fix) = x is a guide to the eye. The error bars indicate sta- 
tistical errors. Systematic uncertainties are 26% and 15% for 
the calculated and measured L^>, respectively. 



the interferometric measurement sequence. 
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